Three-wave mixing of Bogoliubov quasi-particles in a Bose condensate 
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A dressed basis is used to calculate the dynamics of three-wave mixing between Bogoliubov quasi- 
particles in a Bose condensate. Due to the observed oscillations between diflterent momenta modes, 
an energy splitting, analogous to the optical MoUow triplet, appears in the Beliaev damping spectrum 
of the excitations from the oscillating modes. 



Since the experimental realization of Bose-Einstein 
condensation in trapped atomic gases, which serves as 
a mono-energetic and dense atomic source, experiments 
in non-linear atom optics have become feasible. Atomic 
four- wave mixing (4WM) [1], Superradiance [2], and 
matter- wave amplification [3] , [4] are all examples of non- 
linear atom optics, which involve the mixing of several 
atomic fields or the mixing of atomic and electromag- 
netic fields. 

In the case of atomic 4WM, bosonic amplification di- 
rects one of the products of a collision between an ex- 
citation atom and a condensate atom into an, initially, 
largely populated mode [1], [5], [6]. Thus, a new mode 
is macroscopically populated through the second collision 
product. When the energy, e^, of an excitation mode k is 
low compared to the condensate chemical potential, the 
excitation can no longer be described as a free atom mov- 
ing with momentum k, but rather as a collective phonon 
excitation, which involves a large number of atoms [7]. 
Atomic 4WM is therefore inadequate for the description 
of phonon decay. Phonon excitations are described in 
the framework of Bogoliubov theory [8]. As Bogoliubov 
excitations are bosons, bosonic amplification will also di- 
rect one of the damping products of a phonon into an, 
initially, largely populated mode, leading to three wave- 
mixing (3WM) of Bogoliubov quasi-particles. 

In this letter we introduce a basis of states which are 
dressed by the interaction between three, largely pop- 
ulated, modes of Bogoliubov quasi-particles in a Bose- 
Einstein condensate at zero temperature. Using this 
basis of quantum states we calculate the time evolu- 
tion of the system, which exhibits non-linear oscillations 
between the different momentum modes. The use of 
the wave-mixing eigen-basis enables the calculation of 
the system dynamics efficiently and with no approxi- 
mations. In contrast to other theoretical approaches to 
wave-mixing, which lead to non-linear differential equa- 
tions [5] , [9] the dressed state approach turns the problem 
into a linear one, for which propagation in time is triv- 
ial. We show that, due to relative number squeezing, the 
variance in the number difference between the two low 
momentum modes remains constant. 

In analogy to the treatment of spontaneous photon 
scattering in the atom-laser dressed system [10], we treat 



damping into the quasi-continuum of empty modes, dur- 
ing the 3WM dynamics, as transfer between dressed state 
manifolds. Thus, the effects of 3WM on the damping 
process, are calculated. The damping energy spectrum 
is presented. A transition from elastic to inelastic damp- 
ing is observed, which, in analogy to the optical MoUow 
triplet, leads to a splitting of the spectrum into a doublet 
of resonance energies. 

Previously, nonlinear mixing of quasi-particles in a 
Bose condensate was studied, between the two lowest en- 
ergy excitation modes in the discrete regime, where the 
energy of the excited mode is of the order of the energy 
separation between modes [9], [11], [12]. 

Our model system is a homogenous condensate of fi- 
nite volume with Nq atoms in the ground state. 
The Hamiltonian for the system in the Bogoliubov ba- 
sis, taken to the \/Nq order, is given by [13] 
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(1) 



where Hq — ^quNq 



2^:^'^ Ek^fe^^fe + Ek^fc^J^k is the 
part of H which is diagonalized by the Bogoliubov trans- 
formation. A'o is the number of atoms in the k = 
mode. Ek = (S/J + 25?t.) is the recently measured 
Bogoliubov energy [14] and e° is the free particle energy 
e° = h^k^ /2m. g = ATrfi^a/m is the coupling constant, 
a is the s-wave scattering length, m is the mass of the 
atoms, and 6k are the creation and annihilation oper- 
ators respectively, of quasi-particle excitations with mo- 
mentum k. Vk and Uk are the Bogoliubov quasi-particle 
amplitudes. This part of the Hamiltonian describes exci- 
tations in the condensate with momentum k and energy 
£k- 

Hint is the part of the Hamiltonian which is responsible 
for the interaction between excitations 
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Afc,(6+6q6k-q + &+&^q&k). (2) 



The first term in parentheses is referred to as Lan- 
dau damping and is analogous to photonic up conver- 
sion, such as second harmonic generation in optics [15], 
whereas the second term is referred to as Beliaev damp- 
ing, and is analogous to photonic down conversion [16]. 
Akq is the many-body suppression factor 
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Akq = 2Uk [UqUk-q - VqUk-q - UqVk-q) 
-2Vk {VqVk-q — UqVk-q - VqUk-q) ■ 



(3) 



The main damping mechanism from a single, largely 
populated mode k, will be elastic Beliaev damping of 
the excitations into empty modes, which are on a mono- 
energctic surface in momentum space [17]. 

We now consider the case where the condensate is ex- 
cited with N excitations of momentum k and M excita- 
tions of momentum q, jA'^k, Afq, Ok-q), such that k and 
q fulfill the Bragg condition, i.e. Ek = £q + Ck-q- Two- 
photon Bragg transitions can be used to excite the con- 
densate, and populate different momentum modes with a 
variable number of excitations [14]. The amplitude given 
by Hint for Beliaev damping of the k momentum excita- 
tion into two excitations with momenta q and k — q, 

|(iV- l)i^,(Af + l)q,lk-q), will be Vn or ^/M fold 
larger than any of the other damping channels. This 
will result in 3WM into a newly populated k — q mo- 
mentum mode. We assume that iVo >> N,M >> F^o, 
where T is the total Beliaev damping rate of the exci- 
tations into the quasi-continuum of empty modes and 
to is the time of the experiment. Thus, during the ex- 
periment time, 3WM dynamics dominates over damping 
into empty modes. In general the set of N +1 excitation 



Fock states {N — i)^^ , (M -|- i)^ , ik-qy , where i varies 

between and N, spans a degenerate subspace of the 
eigenstates of Hq. Hint couples between pairs of states 
in this subspace which have a difference of 1 in i, and 
can be represented by the (A'^-h 1) x (M + 1) tri-diagonal 
matrix 

{N ~ s,M + s, s\ H,nt \N -i,M + = 



^^,Akq{VN 



- iVM + i + Wi + lSi+i,s + 
y/N -i + WM + iVi5i-i,s)- 



(4) 



When Hint is diagonalized, wc get a new set of TV + 1 
eigenstates, where j varies between 1 and N + 1, 
that are dressed by the interaction. Note that since 
the dressed states arc superpositions of degenerate eigen- 
states of Hq, they are eigenstates of the complete Hamil- 
tonian in Eq. (1). 

The filled circles in Fig. 1 show the calculated energy 
spectrum of the dressed basis for iV = M = 100 in units 
of {g\/N(s/2V)Akq. The degeneracy between the exci- 
tation Fock states is removed by the interaction. The 
spectrum appears to be roughly linear with an average 
energy spacing of dE ~ 2.48\/]V in the above units [18]. 
The hollow circles in Fig. 1 show the value of the energy 
difference between each two dressed states. Due to the 
non-linearity of the problem, energy differences between 
each pair of dressed states are slightly different. 
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FIG. 1. Filled circles show the spectrum of the 
M = N = 100 manifold in units of ^^^Auq. Hollow circles 
show the energy difference between each two dressed states in 
the same units. The spectrum is not linear, and the energy 
spacing varies parabolically around an average of 2.48\/iV. 




FIG. 2. The square of the transfer matrix between the 
excitation Fock state basis and the dressed state basis for 
N = M = IQQ. Darker areas correspond to larger proba- 
bility. The white dashed-dotted line draws the solution to 
\Ej\ = 2^ {N^ — The Fock states are numbered by i. 
The dressed states are numbered by j from lowest to highest 
energy. 

Figure 2 shows the absolute value squared of the trans- 
fer matrix between the bare excitation Fock state basis 
and the dressed basis for A'' = M = 100. Since Hint, af- 
ter being applied a sufficient number of times, can couple 
each two bare states, most dressed states are spanned by 
a large number of bare states [19]. 

In the limit of Af = A'^ ^ oo we can use the tri- 
diagonality of Hint to approximate the dressed state as 
the solution to the differential equation 



9af (x) 



with 



V(l-a;2)a; 
the 



- /o^ aj (x) , 



(5) 



boundary 
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conditions daj {x)/dx = {\j/\/N — aj {x) at a; = 0, 
and dttj {x)/dx = — Xj/^/2N^ aj (a;) at a; = 1. Here 



idx, and aj (x) = {N — i, M + 



dx = 1/iV, a 

is the projection of \j) on the Fock state with i exci- 
tations removed from the k momentum mode. Also, 



A, = EjvN, Ej is the energy eigenvalue of and 

2 - 



fo = 2N\ The main contribution to the dressed state 
superposition of Fock states comes from the two states 



which solve 



i^iJ f 2 

^ Jo 

)x 



0. The dashed-dotted 



line in Fig. 2 draws the solution of this equation. In or- 
der for this equation to have real solutions the resulting 
spectrum must satisfy \Ej\ < 1.24A''^/^. Since there are 
A'' + 1 dressed states we get an average energy difference 
between dressed states of dE ~ 2.48-\/iV, consistent with 
our numerical observation. 




FIG. 3. The expectation value A^'k, k = 0.7 h/^ and 
q = fc/%/2, as a function of time. N = M = 100, for a conden- 
sate of 3 X 10® ®^Rb atoms in the F = 2, mf = 2 ground state. 

The condensate is homogeneous, with a density of 3 x 10^* 
atoms / crn^ . The oscillation frequency is roughly the average 
energy difference in the dressed state spectrum. Insets (a) and 
(b) show Ak and the standard deviation of A^k respectively, 
during a longer time. 

We choose as a model system a condensate of 3 x 
10^ ^"^Rh atoms in the F = 2, m/ = 2 ground state. The 
condensate, which is similar to the experimental param- 
eters of [14] is homogeneous and has a density of 3 x 10^** 
atoms/ cm^. The k = 0.7 H/^ mode and q = mode, 
wIktc ^ is the healing length of the condensate given by 
£, = ^Snna , are populated with a 100 excitations each. 

We start from the above excitation Fock state, written 
as a linear superposition of dressed states, which are ob- 
tained by diagonalizing (4). The state of the system is 
then readily propagated in time by evolving each of the 
dressed states phases according to its energy. Figure 3 
shows iVk, the expectation value of the number of exci- 
tations with momentum k, as a function of time. Excita- 



tions oscillate between the k momentum mode and the q 
and k — q momenta modes. Inset (a) of Fig. 3 shows Nk 
during a much longer time. Since the energy spectrum 
is not precisely linear, at longer times beating between 
the different oscillation frequencies gives rise to a slow 
amplitude modulation of the oscillations. 

Even though we start from an excitation Fock state, 
the system immediately evolves into a superposition of 
Fock states. Inset (b) of Fig. 3 shows the standard devi- 
ation of TVk vs. time. After a time scale which is set by 
the non-linearity of the spectrum, the average value of 
and its standard deviation are of similar size. However, 
the expectation value of the number difference between 
mode q and mode k— q equals M and is constant in time. 
The standard deviation of that difference always remains 
zero, which implies relative number squeezing between 
the two modes. 

Thus far we have discussed the time evolution of the 
system within the dressed state manifold, which is de- 
fined by the initial population of the k and q modes, and 
have neglected damping into empty modes. In analogy to 
the treatment of spontaneous photon scattering as trans- 
fer between dressed states manifolds of the dressed atom- 
laser system [10], we now consider scattering into empty 
modes as transfer between dressed state manifolds. Thus, 
the damping of excitations from mode k is no longer elas- 
tic, but rather carries the energy difference between the 
dressed states among which it occurred. 

Evaluating Hb = ^^iVoEq'^fcg' (^q'^'i^-q'^k) , be- 
tween every pair of states in the two manifolds reveals the 
spectral structure of the damping process. We find that, 
the spectrum of the A^' = N—1, M' = N manifold, which 
has one less energy eigenvalue, is shifted by roughly half 
the energy difference with respect to the N = M energy 
spectrum. Hb significantly couples only dressed states 
with neighboring energies. This results in a structure of 
a doublet in the Beliaev damping spectrum [20]. 

The finite lifetime of the dressed state results from the 
fact that there is a quasi-continuum of N' = N — 1, M' = 
N manifolds, all with an identical energy spectrum, to 
which a dressed state in the N = M manifold can couple 
to via Hb- To determine the width of each transition 
between two dressed states, we use the Fermi golden rule. 
The damping rate between a state ^ in the N, M 



manifold and a state 
is then given by 



N-1,M 



in the N -1, M manifold 



_ 9^ 

h ^ 2y2 



N-1,M 



N,M 



(6) 



5 (£k -I- £o - £q' - £k-q') : 



where eq is the energy difference between ^ and 
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Using momentum conservation the energy con- 
servation (5-function becomes a geometrical condition on 
the angle 6 between k and q' 
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COS (0) = 

1 



(7) 



2kq' 



+ 1 



where momentum is in units of /i/^, and energy is in 
units of gn. Equation (6) is then simplified to 



r = 



j ' 



(£fe + £0 - £q') 



2k 



{ek + eo - £q'f + I 



N-1,M {i\ by, \j)]^^M 

dq'. 



(8) 



The spectrum of each transition is taken as a Lorenzian 
with a width of r/27rA^k and a normalization of T/Nk 
around eq. A^k is the average occupation of mode k for 
the two dressed states involved. Averaging over all of the 
possible transitions between the manifolds, Fig. 4 shows 
the damping spectrum between the N = M = 5 x 10"^ 
and the iV = 5 X 10'"' - 1, M = 5 X 10^ manifolds, for 
the same model system as in Fig. 3 for fc=3.2, 1.6, and 
0.7 h/^, and q = fc/\/2 [21]. A clear doublet structure is 
evident. The inset of Fig. 3 shows the energy-conserving 
surfaces for the two center energies of the k=0.7 h/£, curve 
(solid line), and the energy conserving surface for elastic 
damping from the same mode (dashed line). 
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FIG. 4. Damping spectrum between the N = M = 5 x 10 
manifold and the AT = 5 x 10^ - 1, M = 5 x 10^ manifold, for 
the same condensate as in Fig. 2. The three curves are for 
fc=3.2 (dashed), fe=1.6 (dotted) and fc=0.7 h/^ (solid line), 
q = k/\/2. The inset shows the energy-conserving surfaces 
for the two center frequencies of the k = 0.7 h/^ curve (solid 
line) and the energy conserving surface for elastic damping 
from mode k (dashed line). 

The separation between the peaks in Fig. 4, which 
is equal to the oscillation frequency, decreases with k as 



Akq. The decrease in the width of each peak as a func- 
tion of k comes from two contributions. Firstly, since 
damping into empty modes is incoherent, it decreases as 
jAfeg'l^. Secondly, as k decreases there arc less allowed 
empty modes on the energy-conserving surface. There- 
fore, the doublet structure is more resolved for the lower k 
values. The difference in width between the positive and 
negative energy peaks is mainly due to the difference in 
the number of allowed empty modes on the energy con- 
serving surfaces. For small enough k values, with respect 
to the oscillation frequency, the energy conserving sur- 
face for the negative energy peak completely disappears 
and only the positive energy peak remains. Another con- 
tribution to the width of the resonance in Fig. 4 is due 
to the non-linearity of the dressed state spectrum, which 
gives a slightly different energy for transitions between 
different dressed state pairs. 

Several other mechanisms contribute to the broaden- 
ing of the two peaks which are not included in Fig. 4. 
The fact that only the first scattering event occurs be- 
tween the A'' = M = 5 X 10^ and the AT = 5 x 10^ - 1, 
M = 5 X 10^ manifolds will further broaden the reso- 
nances. Since the energy splitting scales as VlV, in an 
experiment where one scatters dN atoms from mode k, 
this will result in a relative broadening of Accord- 
ing to the same scaling, an initial coherent, rather than 
Fock, state will cause a relative broadening of The 
condensate finite size or inhomogcncous density profile, 
will further contribute to the width of the resonance. We 
estimate that for the experimental parameters of [14] a 
doublet structure in the Bcliacv damping spectrum can 
be resolved. Experimentally, the energy doublet can be 
observed by computerized tomography analysis of time 
of flight absorption images of the 3WM system [22]. 

Another process which will transfer the system be- 
tween manifolds is that of a two-photon Bragg transition 
from the condensate to cither the k or q modes [14] . The 
line-shape of the Bragg transition into the k mode also 
splits into a doublet structure [23]. 

The dressed state formalism can be applied to other 
processes which involve bosonic amplification. We 
have applied this method to a strongly Bragg-driven 
condensate, with large momentum transfer, by diago- 
nalizing Hint = ^^^^ {a^ao + a^) , is the two- 
photon Rabi frequency, and to atomic 4WM with 
large relative momentum by diagonalizing Hint = 

2y (^Q^k '^d'^k-qaq + ftk-qflq '^kO'O^ • Both systcms show 
a doublet structure in the spectrum of collisions from a 
certain momentum mode, due to rapid oscillations in the 
population of that mode [23]. 

In conclusion, we calculate the wave-mixing dynam- 
ics between three, low k, Bogoliubov quasi-particles in a 
Bose condensate. The Hamiltonian of this system is di- 
agonalized to the next order in The resulting basis 

of dressed states allows for the efficient, linear, propaga- 
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tion of the system in time. Non-linear oscillations be- 
tween the different momentum modes are observed. Rel- 
ative number squeezing between the q and the k — q 
momentum modes is shown. Bcliacv damping of excita- 
tions from these modes is treated as a transfer between 
dressed states manifolds. The damping process is shown 
to become inelastic and, similarly to the optical MoUow 
triplet, exhibits a doublet-like energy spectrum, which is 
more resolved for the low k values. 
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5 X 10^, we diagonalizc smaller N.M matrices and scale 

the resulting energies by a factor of \/iV. 
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